QCD: Restoration of Chiral Symmetry and Deconfinement for Large Nj 
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Exploiting the recent lattice results for the infrared gluon propagator with light dynamical quarks, 
we solve the gap equation for the quark propagator. We thus extract the order parameters for chiral 
symmetry breaking (constituent quark mass, quark-antiquark condensate and the pion leptonic de- 
cay constant) as well as confinement (complex conjugate mass-like singularities) which is intimately 
tied with the analytic properties of the Schwinger functions of quantum chromodynamics (QCD). 
We demonstrate that there are clear signals of chiral symmetry restoration and deconfinement when 
the number of light quark flavors exceeds a critical value of Nf ~ 8. 
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Quantum chromodynamics (QCD) with large number 
of massless fermion flavors has seen a resurgence of inter- 
est due to its connection with technicolor models, origi- 
nally proposed by Weinberg and Susskind pLj, which fall 
into the category of Beyond the Standard Model Theo- 
ries. They possess intrinsically attractive features. They 
do not resort to fundamental scalars to reconcile local 
gauge symmetry with massive mediators of interaction 
and have close resemblance with well-studied fundamen- 
tal strong interactions, i.e., QCD. However, their simple 
versions do not live up to the experimental electroweak 
precision constraints, in particular the ones related to 
flavor changing neutral currents. Walking models con- 
taining a conformal window and an infrared fixed point 
can possibly cure this defect and become phenomenologi- 
cally viable j5] . This scenario motivates the investigation 
of QCD for similar characteristics. One looks for such be- 
havior of QCD for large number of light flavors albeit less 
than the critical value where asymptotic freedom sets in, 
i.e., N^'^ = 16.5, a Nobel prize winning result known since 
the advent of QCD, Just as Nf dictates the peculiar 
behavior of QCD in the ultraviolet, we expect it to de- 
termine the onslaught of its emerging phenomena in the 
infrared, i.e., chiral symmetry breaking and conflnement. 

Whereas the self interaction of gluons provides anti- 
screening, the production of virtual quark-antiquark 
pairs screens and debilitates the strength of this interac- 
tion of non abelian origin. For real QCD, light flavors are 
small in number and hence yield to the gluonic influence 
which triggers conflnement and chiral symmetry break- 
ing. One needs to establish if there is another critical 
value Nj^ < N^^ which can sufficiently dilute the gluon- 
gluon interactions to restore chiral symmetry and decon- 
flne color degrees of freedom. Such a phase transition 
lies at the non perturbative boundary of the interactions 
under scrutiny and hence we cannot expect to extract 
sufficiently reliable information from multiloop calcula- 
tions of the QCD /^-function. Purely non perturbative 
techniques are required to tackle the problem. Lattice 
studies in the infrared indicate that just below N^^ , chi- 
ral symmetry remains unbroken and color degrees of free- 



dom are unconfined [4 . Below this conformal window, 
for an 8 < Nj^ < 12, the evolution of the beta func- 
tion in the infrared is such that QCD enters the phase of 
dynamical mass generation as well as confinement. 

In continuum, Schwinger-Dyson equations (SDEs) of 
QCD provide an ideal framework to study its infrared 
properties, [5_. These are the fundamental equations of 
any quantum field theory, linking all its defining Green 
functions to each other through intricately coupled non- 
linear integral equations. As their formal derivation 
through variational principle makes no appeal to the 
weakness of the interaction strength, they naturally con- 
nect the perturbative ultraviolet physics with its emerg- 
ing non perturbative properties in the infrared sector 
within the same framework. The simplest two-point 
quark propagator is a basic object to analyze dynami- 
cal chiral symmetry breaking and confinement. Within 
the formalism of the SDEs, the inverse quark propagator 
can be expressed as 

S-\p) = Z2(n-P + ™) + S(p) , (1) 

where S(p) is the quark self energy: 

^^P^^^^ J (0iff'A^.b-9)y7M^(9)r^(9,p), (2) 

where Zi — Zi(^^,A^) and Z2 — Z2(^^,A^) are the 
renormalization constants associated respectively with 
the quark-gluon vertex and the quark propagator. A 
is the ultraviolet regulator and fi is the renormalization 
point. The solution to this equation is 

where Z{p^,y?) is the quark wave function renormal- 
ization and the quark mass function M(p'^) is renor- 
malization group invariant. This equation involves the 
quark-gluon vertex r°(g,p) and the gluon propagator 
A^^(p- g). 

As a consequence of a patient effort, spanning several 
decades to unravel gluon propagator A^j^ in the infrared. 
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lattice as well as SDE studies have finally converged on 
its massive or so called decoupling solution, see for ex- 
ample [B]. After the gluon propagator solution in the 
quenched approximation has been chiselled, we now have 
the first quantitatively reliable glimpses of its quark fia- 
vor dependence by incorporating |29] Nf = 0, 2, 4 dy- 
namical quark flavors in the QCD action. In Ref. |7j, the 
authors studied the gluon propagator with light quark 
masses ranging from 20 to 50 [MeV], while the strange 
quark is roughly set to 95 [MeV] and the heavy charm 
to 1.51 [GeV] (in MS at ^ 4 [GeV^]). Very recently 
also, the authors of Ref. [S] developped a "partially un- 
quenched" approach to incorporate flavor effects in the 
gluon SDE and yielded results in agreement with those 
in Ref. [7 . In any case, this two-point function serves 
as a crucial input to study the quark propagator. The 
only other ingredient is the three-point quark-gluon ver- 
tex r"(g,p). Signiflcant advances have been made in pin- 
ning it down through its key attributes in the ultraviolet 
and infrared domains [S]. More recently, the seeds of 
the most general ansatz for the fermion-boson vertex ap- 
peared in [10_ and its full blown extension was presented 
in |11| . Signiflcantly, this ansatz contains nontrivial fac- 
tors associated with those tensors whose appearance is 
expressly driven by dynamical chiral symmetry breaking 
in a perturbatively massless theory. This novel feature 
enables a direct and positive comparison with the best 
available symmetry preserving solutions of the inhonio- 
geneous Bethe-Salpeter equation for the vector vertex. 
The encouraging outcome indicates that this model is 
likely to provide a much-needed tool for use in Poincare- 
covariant symmetry-preserving studies of hadron electro- 
magnetic form factors. Furthermore, given the general 
nature of constraints and the simplicity of the construc- 
tion, a straightforward extension of this approach is ex- 
pected to yield an ansatz adequate to the task of rep- 
resenting the dressed-quark-gluon vertex. Before this is 
achieved, we restrict ourselves to an effective though effi- 
cacious approach. Following the lead of Maris et. al. |12| . 
we employ the following suitable ansatz which has suf- 
ficient integrated strength in the infrared to achieve dy- 
namical mass generation: 

Z^g'A,,{p^q)r,{p,q)^gl^ A^,(p~q)^7, , (4) 



where 



KM 



D{q^ 



(5) 



The effective coupling g^s is chosen to correctly repro- 
duce the static as well as dynamic properties of mesons 
below 1 GeV and reproduce perturbation theory in the 
ultraviolet, see for example review [13] and references 
therein. Moreover, our modern understanding of the 
flavor dependence of the gluon propagator provides us 
with the solid basis to use the following non perturbative 



model [Ti) : 



z(m') q^jq^ + M^) 

q^ + q^ (m2- 13.92 (A2)/24) 



(6) 







to describe the gluon dressing renormalized in MOM 
scheme at q^ = fj?. This model is based on the tree- 
level gluon propagator obtained with the Renormalized 
Gribov-Zwanziger (RGZ) action [TS] which have been 
shown to describe properly the lattice data in the infrared 
sector (see refs. [HI IT7]V The overall factor ^(/i^) is 
introduced to guarantee the multiplicative MOM renor- 
malization prescription, namely, D{^^) = 1, and implies 
no physical consequence as the effective coupling, ^cff , is 
furtherly adjusted to reproduce properly the meson phe- 
nomenology. The mass parameter from Eq. (|6| is 
related to the condensate of auxiliary fields only emerg- 
ing to preserve locality for the RGZ action, i.e., when 
incorporating the horizon condition to the action. {A"^) 
is the dimension-two gluon condensate that has been 
elsewhere very much studied [16]). niQ is another mass 
scale, m§ = z{fM^) \imq2^Q q'^ /D{q^). It is related to 
the above mentioned condensates and to the Gribov- 
Zwanziger parameter determined by the well-known hori- 
zon condition. To determine its mass parameters, Eq. ^ 
can be fitted to the gluon propagator lattice data ana- 
lyzed in Ref. [7] with the goal of scrutinizing the quark 
flavor effects on ghost and gluon propagators. These 
data had been orginally obtained from gauge conflgura- 
tions simulated at several lattices with Nf=2 [18j and 
Nf=4 [19_ mass- twisted lattice flavors [20] (generated 
within the framework of ETM collaboration [^,'5^), and 
with Nf = ||23j. As can be seen from Fig. [l] we thus 




FIG. 1: Parameters < A'^ > and l/mg in terms of the 
numbers of flavors and the linear fits yielding the results of 
Eq. |7|. 

obtain the following behavior in terms of the number of 
fermion flavors for the mass parameters: 



mo = 1.011 (9.161 - 
g'^{A^) = 0.474 (16.406 
m2 = 4.85 [GeV^] ; 



[GeV] , 



Nf) [GeV^] 



(7) 
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where is required not to depend on the number of 
fermion flavors. This is a reasonable approximation, 
fairly confirmed by a free fit. These parameters, plugged 
into Eq. (|6|, result in the giuon propagator prediction, for 
any fiavor number, shown in Fig. [2] to describe properly 
the lattice data for Nf —0,2,4 of Ref. [7J. Thus, we can 
efficaciously model the dilution of the gluon-gluon inter- 
actions with increasing flavor number in order to study 
the chiral restoration mechanism. 
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FIG. 2: Gluon propagator in terms of momenta, as given 
from Eqs.( |6|7| , for different number of fermion flavors. The 
points correspond to the same "unquenched" (A^/=2,4), and 
quenched (A'^/=0) lattice data analyzed in Ref. [7]. 

We can now employ the gap equation to provide quan- 
titative details of chiral symmetry breaking in terms of 
the quark mass function for an increasing number of light 
quarks. We take the effective coupling ^cff to be inde- 
pendent of Nf which is justifled by the results of [7J (see 
Eq.(5.2)) which suggest that an effective coupling can be 
constructed such that there is an absence of any flavor 
dependence in the infrared region, more precisely start- 
ing from < 1 GeV^. Note that we have not considered 
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FIG. 3: The quark mass function diminishes height for in- 
creasing light quark flavors. Above Nf = 8.33685, only the 
chirally symmetric Wigner solution exists. 

the flavor dependence which would arise from the quark- 
gluon vertex. No explicit handle on this dependence is 



available at the moment. Within the abelian theory of 
QED, restrictions imposed by the all order multiplicative 
renormalizability of the photon propagator may provide 
a handle on the transverse part of the electron-photon 
vertex, see the last reference in |9j. A consequent non 
perturbative construction of such a vertex with imprints 
of the massless charged fermion flavors and its subse- 
quent extension to QCD is still not available. Once the 
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FIG. 4: Quark pole mass in the Euclidean space clearly 
demonstrates that chiral symmetry is restored above a criti- 
cal number of quark flavors. The solid line is the mean-field 
scaling fit described by Eq. Q. 

quark mass function is available for varying light quark 
flavors, one can investigate any of the interrelated order 
parameters, namely, the Euclidean pole mass deflned as 
™dyn + M^{p'^ = TO^yjj) = 0, the quark-autiquark con- 
densate which is obtained from the trace of the quark 
propagator or the pion leptonic decay constant de- 
flned through the Pagel-Stokar equation [21] , or through 
considering the residue at the pion pole of the meson 
propagator. Each of these quantities involves the quark 
wave- function renormalization, the mass function and/or 
its derivatives and is hence calculable from the solution 
for the full quark propagator. Moreover, these order pa- 
rameters can help locate the critical number of flavors 
above which chiral symmetry is restored. We choose to 
present the Euclidean pole mass of the quark, shown in 
Fig. |4] At a critical value of about N^ « 8, chiral sym- 
metry breaking ceases to exist. The phase transition ap- 
pears second order, described by the following mean fleld 
behavior (shown by the solid line in Fig. |4]) : 

mdyn ~ ^JN^/ ~ Nf . (8) 

It has been established that conflnement is related to the 
analytic properties of QCD Schwinger functions which 
are the Euclidean space Green functions, namely, propa- 
gators and vertices. One deduces from the reconstruction 
theorem |25j that the only Schwinger functions which can 
be associated with expectation values in the Hilbert space 
of observables; namely, the set of measurable expectation 
values, are those that satisfy the axiom of reflection posi- 
tivity. When that happens, the real-axis mass-pole splits. 
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FIG. 5: Spatially averaged Euclidean space 2-point Schwinger 
function A(t) develops oscillations for large times which cor- 
responds to the non-existence of asymptotically stable free 
quark states. For sufficiently large values of Nf, the first 
minimum of these oscillations is pushed all the way to infin- 
ity, thus ensuring the existence of a pole on the time-like axis, 
a property of free particle propagators. 



A{t) 



(2^ 



^j(p4t + p-x) 



(9) 



to study the analytic properties of the quark propaga- 
tor; where crs(p^) is the scalar term for the quark prop- 
agator in Eq. ([3|, that can be written in terms of the 
quark wavefunction renormahzation and mass function 
as Z{p'^,fi'^)M{p'^)/{p'^ + M{p'^)). One can show that if 
there is a stable asymptotic state associated with this 
propagator, with a mass to, then A(f) ^ e~™*, whereas 
two complex conjugate mass- like singularities, with com- 
plex masses ^ ^ a±ib lead to an oscillating behavior of 
the sort A(t) ^ e-"*cos(6t + 5) for large t, Fig. [s] 

analyzes this function for varying Nf. The existence of 
oscillations clearly demonstrates that the quarks corre- 
spond to a confined excitation for small Nf. With in- 
creasing Nf, the onslaught of oscillations moves towards 
higher values of t and eventually never takes place above 
a critical Nf when quarks deconfine and correspond to a 
stable asymptotic state. 



moving into pairs of complex conjugate singularities. No 
mass-shell can be associated with a particle whose prop- 
agator exhibits such singularity structure. We define the 
following Schwinger function: 




As an order parameter of confinement, we therefore 
employ v{Nf) = l/Ti{Nf), where Ti{Nf) is the location 
of the first zero, j57j. This order parameter vanishes when 
confinement is lost. It is notable that when the dynami- 
cally generated mass approaches zero, i'{Nf) diminishes 
rapidly, as can be seen in Fig. |6] This highlights the in- 
timate connection between chiral symmetry restoration 
and deconfinement. In fact, within the numerical accu- 
racy of our work. Fig. [6] shows Nj to be the same for 
both the transitions. 



Therefore, the SDK analysis of the latest lattice results 
for the quark fiavor dependence of the gluon propagator 
in the infrared hint towards chiral symmetry restoration 
and deconfinement in QCD when the number of light 
quark fiavors exceeds a critical value of N'j'^ ss 8. Notice- 
ably, roughly the same number is also an indication of 
modern lattice QCD analysis |28] . 



FIG. 6: The order parameter for confinement i^{Nf) = 
l/ri{Nf), where Ti{Nf) is the location of the first zero 
of Eq. |9|. Comparison with Fig. [4] suggests that quarks get 
deconfined when chiral symmetry is restored. 
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